SU(4) Kondo Effect in Carbon Nanotubes 
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We investigate theoretically the non-equilibrium transport properties of carbon nanotube quantum 
dots. Owing to the two-dimensional band structure of graphene, a double orbital degeneracy plays 
the role of a pseudo-spin, which is entangled with the spin. Quantum fluctuations between these 
four degrees of freedom result in an SU(4) Kondo effect at low temperatures. This exotic Kondo 
effect manifests as a four-peak splitting in the non-linear conductance when an axial magnetic field 
is applied. 
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Introduction. — Carbon nanotubes (CNTs) exhibit a 
good deal of remarkable transport phenomena including 
quantum interference 0, Luttinger liquid features or 
spin polarized transport Q . Finite- length CNTs behave 
like quantum dots (QDs) and thus show Coulomb block- 
ade jj, |]| and Kondo physics . Interestingly, the rich- 
ness of the band structure of CNTs and the feasibility to 
attach new materials as electrodes, e.g., ferromagnetic |3| 
or superconducting contacts pj , allows us to explore new 
aspects of the Kondo effect, one of the most fundamental 
topics in condensed matter physics. 

The electronic states of a CNT form one-dimensional 
electron and hole sub-bands. They originate from the 
quantization of the electron wavenumber perpendicular 
to the nanotube axis, k±, which arises when graphene 
is wrapped into a cylinder to create a CNT. By symme- 
try, for a given sub-band at fej_ = ko there is a second 
degenerate sub-band at k± = — ko. Scmiclassically, this 
degeneracy corresponds to the clockwise (O) or counter- 
clockwise (O) symmetry of the wrapping modes. Linear 
transport measurements in the Coulomb blockade regime 
reveal a distinct four-fold shell filling pattern owing to the 
orbital and the spin degeneracies [5|. 

In this Letter we combine several theoretical ap- 
proaches, scaling theory, numerical renormalization 
group (NRG) non-crossin g ap proximation (NCA) [9|, 
equation-of- motion (EOM) |l(J methods, to present 
a unified picture of low-temperature, non-equilibrium 
transport through CNT quantum dots in the presence 
of magnetic fields. We show that quantum fluctuations 
between the four states {Of, 0|, Of, 01} may dominate 
transport at low temperatures provided that both the 
orbital and spin indexes are conserved during tunneling. 
This leads to a highly symmetric SU(4) Kondo effect, and 
hence an enhanced Kondo temperature, in which the spin 
and the orbital degrees of freedom are totally entangled 
fl"i| . We also point out that the orbital degeneracy in the 
dot itself is not enough for having SU(4) Kondo physics. 
In general, SU(2) Kondo physics is possible. We show 
that neither an enhanced Kondo temperature nor linear 
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FIG. 1: A schematic of couplings between the degenerate 
orbitals [m = 1(2)] with conduction channels [L(R),m]. 



conductance measurements can distinguish between the 
two effects. Instead, the non-linear conductance in the 
presence of a parallel magnetic field shows a four-peak 
structure, with different splittings for the spin and the or- 
bital sectors, which unambiguously signals SU(4) Kondo 
physics. Our theoretical results are in good agreement 
with recent experiments by Jarillo-Herrero et al |l2^ . 

Model. — We study a quantum dot CNT coupled to 
left (L) and right (R) electrodes. We consider the case 
where the QD has two (nearly) degenerate localized or- 
bitals [labeled by the quantum number m = 1,2 where 
1(2) denotes O (0) orbital mode, respectively]. The pres- 
ence of an axial magnetic field (-Bn) lifts both the orbital 
and spin degeneracies: A parallel magnetic field induces 
an Aharonov-Bohm phase 27r*E>/ *E»o, where $ = ndf/AB^ 
is the flux threading the nanotube, $o = h/e is the 
flux quantum and dt the tube diameter. This Aharonov- 
Bohm flux shifts the allowed k± and the orbital degener- 
acy is lifted by an amount ±A or b = iedtVpBu/A, where 
Vf if the Fermi velocity The states near the en- 

ergy gap correspond to semiclassical orbits which have 
an orbital magnetic moment /i or b = ed t vp/A. Thus, 
an axial magnetic field leads to an energy shift A or b = 
Morb^||- The states further split due to the Zeeman en- 
ergy ±Az/2 = ±<?/iB-B|i/2, where fiB is the Bohr mag- 
neton and g ~ 2 is the g-factor in nanotubes. /i or b 
scales with the CNT diameter and is typically one or- 
der of magnitude larger than the Bohr magneton [lil FLU ] 
resulting in a stronger influence of Bu on the orbital 
sector. Thus, the single particle energy e ma associ- 
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ated with the orbital m and the spin a is given by: 
f-mo- = (d + A orb (6 m ,i - 8 m<2 ) + (A z /2)(S a ^ - 8^1). The 
QD is then described by the Hamiltonian 

ro=i,2<x=T,| 

where n = J2 m a d m a-dm,a is the occupation and U is the 
Hubbard-like on-site interaction (we focus on the regime 
where the QD is occupied by a single electron) . The two 
leads are modelled as: 

-^C = E 22 ea > kC <x,k,m,<T C o<,k,m,cT , (2) 

a£L,Rm=l,2 k,a 

Without loss of generality, we assume that there are two 
distinguished (groups of) channels m = 1, 2 in each lead. 
The coupling between the leads and the dot is described 
by the tunneling Hamiltonian of the following form 

*=EE V Z,™> (< fc , m , CT <W + h.c) , (3) 
a,fc,<7 m,m' 

where the tunneling amplitudes read V 7 " m , = [Vb^m.m' + 
Vx(l — <5m,m')]/v / 2, (for simplicity, we ignore the k- and 
CT-dependence of the tunneling amplitudes). In the most 
general case, they describe two types of tunneling pro- 
cesses (see Fig.^): (i) those in which the orbital quantum 
number is conserved, denoted by V^ x = V 2 a 2 = Vo/y/2 
and (ii) those events accounting for mixing (cross cou- 
pling) y ? " m , — Vx / V2 with m 7^ m! . To gain more phys- 
ical intuition when Vx we rewrite the total Hamilto- 
nian Jj? = J4?d + Jf?c + Mr in terms of symmetric (even) 
and antisymmetric (odd) combinations of the orbital 
channels. Before, we simplify the algebra by perform- 
ing a canonical transformation, cum k m & — (a,h,m,a ± 
bk,m,a)/V^, such that the resulting Hamiltonian contains 
only a single lead, )/v / 2, 
with two channels m = 1,2. Next, we apply the even- 
odd transformation a^mVer = (ck c a ± ick a)/V% and 
di(2)a = {dea ± id oa )/y/2, such that Jf? reads, 

= Y e ku c k Ut<T Ck u ,<T+ Y € v?di,<jdv,<j+U{n-n g ) 

+ ye E ( C U<*~ + h.c) +V°Y (4 ot ,do, + h.c) , 

(4) 

with V e = V + V x and V° = V - V x (note that e kv 
and e„ CT remain invariant under this transformation). In 
the absence of cross coupling terms Vx = both, the 
even and odd orbitals are equally coupled to the QD 
(V e = V° = Vq). Thus, at small energies the effective 
model leads to SU(4) Kondo physics, see below. However, 
for the maximal mixing, i.e, Vx = Vq, the even orbital 
is doubly coupled to the dot (V e — 2Vq) whereas the 



odd orbital becomes uncoupled (V° = 0). Here, SU(2) 
Kondo physics arises owing to spin fluctuations in the 
even orbital channel. 

Effective Kondo model. — Let us now substantiate our 
previous arguments by examining the low-energy prop- 
erties of the system. After performing a Schrieffer-Wolf 
transformation to Jif the effective Hamiltonian reads: 

J? K = Jf C + -j [S ■ + S • {^(TT Z ljj)T x ] 

+ [S ■ (V-Vt^V) ■ T- 1 + (^r-V) ■ T x ] 
+ ■h^ T ^)T x + J -± [S • (^(TT z ip) + S • (^ct^)T x \ 
- J 5 T X + n B gB • S + /i orb B • T, (5) 

where $\ = [c\ kA c\ M c\ M 4 >fe J are the field opera- 
tors for the leads in the even-odd basis and ip = ^ fc ipk- 
The field operator for the QD i>\ = [d\ ^ d\ , d' + d\ J 
defines the orbital pseudo-spin and the spin operators 
given by T = ip^Ttpd and S = ip^atpd respectively, with 
t and a being the Pauli matrices in the orbital pseudo- 
spin and spin spaces, respectively. Here, T^(t^) de- 
notes T y + T z (t v + t 2 ). Using renormalization group 
(RG) arguments, the effective coupling constants in the 
[/-large limit are given initially by J\ = J3 = ,jVV 2 /e^, 
■h = Ji{\V \ 2 - \V X \ 2 )/V 2 , J 4 = ,h\2V V x \/V 2 , where 
V 2 = \V \ 2 + \Vx\ 2 {yy is the degeneracy). J 5 is not renor- 
malized in the RG procedure and does not flow into the 
strong coupling regime. At zero magnetic field with only 
spin and orbital conserving tunneling processes (Vx = 0), 
J 1 = J2 = J3 = J while J4 = J5 = 0. The corresponding 
Hamiltonian is reduced to the SU( t /( / =4) Kondo model 
where the spin S and the orbital pseudo-spin T are en- 
tangled (last term of the equation), 

^su(4) = + (J/4)[S • tyW) + WW) ■ T 

+ S • (ip^a-Tip) ■ T]. (6) 

The scaling equations are reduced to a single equation 

dJ/InD = -jVp J 2 , (7) 

where po is the density of states (DOS) in the leads and D 
is the bandwidth, resulting in the exponentially enhanced 
Kondo temperature w Dexp[l/(.yVp J)] com- 

pared with T^ U(2) of the single-level SU(^( / =2) model. 

In the other limiting case (Vx = Vq), the correspond- 
ing Kondo-like Hamiltonian [Eq. J5J at B — with J\ — 
J3 = J4 = 2Jf I Vb| 2 /erf and J2 — 0)] involves only spin 
fluctuations in the doubly degenerate (,jV=2) even or- 
bital. It gives rise to what we call a two-level (TL) SU(2) 
model: ^tlsu(2) = + JS e • (^ e a^ e ) (1 + T x ) + 
(J/A)(iljl^ e )T x - J 5 T X . The RG equation for J is 

dJ/lnD = 2,yVp J 2 . (8) 
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This model corresponds to an SU(2) Fermi liquid [l5| 
with a Kondo temperature yTLSUC 2 ) wn j cn j s ^ e 



same 

as T^ V(4> (compare Eqs. (0 and ©, the factor 2 comes 
from the doubling of the coupling 2Vq (lfijl- 

The scaling arguments above are confirmed by the 
NRG studies of the spectral density A ma (u)) for the lo- 
calized level ma. At B\\ = 0, the spectral density shows 
a peak near the Fermi energy, corresponding to the for- 
mation of the SU(4) Kondo state; see Fig. [S] (solid line 
in the right inset). The peak width, which is much 
broader than that for the SU(2) Kondo model (dotted 
line), demonstrates the exponential enhancement of the 
Kondo temperature mentioned above. Another remark- 
able effect is that the SU(4) Kondo peak shifts away from 
w = Ep = and is pinned at lu w T^^ 4 \ This can be 
understood from the Friedel sum rule [l?J which, in this 
case, gives 8 — it /A for the scattering phase shift at Ep. 
Accordingly, the linear conductance at zero temperature 
is given by Go = 4(e 2 /ft.) sin 2 8 = 2e 2 jh. It is interest- 
ing to recall that the Friedel sum rule gives the same 
linear conductance also for the TL SU(2) Kondo model. 
Thus, neither the enhancement of the Kondo tempera- 
ture nor the linear conductance, can distinguish between 
the SU(4) and the TL SU(2) Kondo effects. This can 
only be achieved by studying the influence of a parallel 
magnetic field, which we do now. 

SU(4) Kondo model at finite field. — Because of the un- 
derlying SU(4) symmetry, the orbital pseudo-spin should 
behave the same way as the real spin. In particular, the 
lift of the pseudo-spin degeneracy will split the Kondo 
peak (as long as the lift is larger than the Kondo tem- 
perature) just like the Zeeman splitting of the real spin 
does. The only difference is that pseudo-spin is more 
susceptible to the magnetic field than the real spin since 
Morb ^> \ib (see above). Therefore, at sufficiently large 
fields (2A orb > A z > T^ U(4) ), one has four split-Kondo 
peaks at u ±2A or b and u> « ±Az; see Fig. (left 
inset). At moderate fields such that 2A or b > » 
T^- > Az, one can have a three-peak structure; see 
Fig- ID (dashed line). The lifted degeneracy in the orbital 
pseudo-spin gives two side-peaks at u> w ±2A or b while 
the spin still retains a Kondo effect and gives the central 
peak. The central peak (which is now at w = 0) corre- 
sponds to a conventional SU(2) Kondo effect and hence 
is much narrower than the central resonance for Bu — 0. 

The above features of A to t(uj) at equilibrium 
are directly reflected in the non-linear conductance, 
Q = dl/dV, an experimentally measurable quan- 
tity. The current through the system I can be 
expressed in terms of the local DOS ^^| : I = 
(e/*)E„=u £«=!,:» /dw - f R (uj)}[T (uj) + 

8T(uj)}A ma (uj), where r (u/) = np (us)\Vo\ 2 , 8r{uj) = 
irpo(uj)\Vx | 2 , and /l(r) is the Fermi function for the 
left (right) lead. The NRG procedure is valid only at equi- 
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FIG. 2: (Color online) NRG results for the total spectral den- 
sity A to t(u}). The solid line is for the case of A or b = Az = 0, 



the dashed line for A 
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A z w 0, and the dash-dotted line for A orb = 4A Z = ST t 
(Parameters: e d = -10r , U = 200r , To = 0.01D," and 
ST = 0). Left inset: Zoom of the four peak splitting in 
-Atot(w)- Right inset: Comparison of the SU(4) (solid line) 
and the single- level SU(2) (dotted line) Kondo models. (Pa- 
rameters: e d = ~5L , U = 500r , r = 0.02D, and 5F = 0). 



librium so we need a method capable to produce the non- 
equilibrium DOS and the non-linear current. We choose 
to use a combination of the NCA and the EOM meth- 
ods Figures Of a) and (b) display the local DOS and 
the differential conductance, respectively, for several B» . 



We take A or b ranging from A or b = 0.5T 
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to A orb = 

1.5T^ U ^ with /i or b = fO^s leading to Az = A or b/5. For 
the lowest magnetic field, A to t{uj) exhibits two splittings, 
namely the orbital and the Zeeman splittings with peaks 
atwsi ±2A or b and to f=a ±A^, respectively, in excellent 
agreement with the NRG calculations. As A or b enhances 
beyond, a substructure arises in the outer peaks. These 
new side-peaks at ui ~ ±2A or b ± Az correspond to the 
simultaneous spin- flip inter- orbital transitions. Increas- 
ing further B\\ these side-peaks in the DOS get better 
resolved 0|. All these features are also present in the 
differential conductance plotted in Fig.OJb). 

Two-level SU(2) Kondo model at finite field. — As al- 
ready shown, both the SU(4) and TL SU(2) Kondo mod- 
els lead to similar results for the Kondo temperature 
and the linear conductance. Nevertheless, the inherent 
physics is completely different. In the highly symmet- 
ric SU(4) Kondo model, the orbital pseudo-spin and the 
real spin are indistinguishable (and screened simultane- 
ously at zero field). On the contrary, in less symmetric 
multiple-level SU(2) Kondo models, the spin should be 
clearly distinguished from the orbital sector (tunneling 
processes preserve only the spin). This becomes clear at 
finite magnetic fields: When B« lifts the orbital degener- 
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FIG. 4: (Color online). NRG results for the total spectral 
density comparing (a) the SU(4) Kondo model (ST = 0) and 
(b) the two-level SU(2) Kondo model (ST = To)- Parameters 

are A orb = 5A Z = 15T£ U(4) (« 0.2r ), e d = -10r , U = 
1000r , and T = 0.0LD. 



FIG. 3: (Color online). NCA+EOM results for the SU(4) 
Kondo model, (a) .Atot(w) and (b) dl/dV versus eV for dif- 
ferent magnetic fields ranging from A orb = 0.5T^ U(4) (bot- 
tom) to A or b = 1.5TjL u ^ (top). (The curves are shifted 
vertically for clarity). When Bu ^ 0, the Kondo resonance 
splits due to the removal of both spin and orbital degenera- 
cies. The rest of parameters are: = — 4Fo and T = 0.003Fo 
with ST = 0. The quantum dot is symmetrically coupled to 
two leads (r_L = Tr = To) consisting of Lorentzian bands 
of width 2D — 20ro. Inset: schematic representation of the 
allowed transitions. 1) intra-orbital with spin- flip, 2) inter- 
orbital without spin-flip and 3) inter-orbital with spin-flip. 



acy, only the lower orbital level is occupied so the physics 
is essentially that of a sing lc level Kondo model 161. The 
Kondo resonance peak gets narrower with increasing Bu . 
Since B\\ also breaks the spin degeneracy, the resulting 
DOS displays the usual Zeeman splitting. Overall, there 
are only two peaks around Ep (Fig- EJ ■ 

Summary. — We have demonstrated that quantum 
fluctuations between the orbital and the spin degrees 
of freedom in carbon nanotube QDs result in an SU(4) 
Kondo effect at low temperatures. This exotic Kondo ef- 
fect manifests as a four-peak splitting in the non-linear 
conductance when an axial magnetic field is applied. 
Similar effects may appear in other systems like verti- 
cal dots |2fj| where the orbital quantum number is pre- 
served during tunneling. Recent transport experiments 
in CNTs 01 dearly support our theoretical findings. 
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